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ABSTRACT 

A  septum  and  an  Impedance  matching  post  are  used  as  a  beam  equalizer 
in  an  open-ended  waveguide- feed  for  reflectors  used  in  satellite  communica¬ 
tions  systems*  The  performance  of  this  design  over  a  frequency  band  is 
evaluated  using  a  spectral  domain  approach*  The  computed  radiation 
patterns  in  the  E-  and  H-planes,  as  veil  as  the  results  for  the  impedance 
match,  ire  presented  in  the  paper* 
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I.  Introduction 


Rectangular  waveguide  array  feeds  for  reflector  antennas  play  an 
important  role  in  the  design  of  satellite  communication  systems.  To  make 
the  radiation  pattern  more  symmetric  in  the  E-  and  H-planes  of  the  feed, 
a  beam  equalizer  is  needed.  The  design  used  in  this  case  is  a  septum  placed 
across  the  mouth  of  the  waveguide  such  that  the  aperture  distribution  is 
reshaped  to  satisfy  the  new  boundary  conditions  imposed  by  the  septum. 
Consequently,  the  H-plane  radiation  pattern  is  narrowed  to  approach  the 
E-plane  pattern,  thereby  achieving  the  beam  equalizing  effect.  However, 
the  introduction  of  such  a  septum  creates  an  impedance  mismatch  problem  for 
the  feed*  To  alleviate  this  problem,  a  matching  post  is  placed  behind  the 
septum  so  that  the  reflection  back  into  the  waveguide  is  minimized. 

The  performance  of  this  design  over  the  desired  frequency  band  is 
evaluated  using  a  spectral  domain  approach,  or  more  specifically,  Galerkin's 
method  applied  in  the  spectral  domain  [1].  The  scattered  fields  on  both 
sides  of  the  beam  equalizer  are  represented  in  terms  of  their  Fourier 
transforms  or  spectra  which  can  be  related  to  the  induced  surface  currents 
on  the  septum  and  the  post.  These  unknown  induced  currents  are  expanded  in 
terms  of  known  basis  functions  and  unknown  coefficients.  A  matrix  equation 
for  the  unknown  coefficients  is  derived  by  applying  the  boundary  conditions, 
and  the  moment  method  is  then  employed  in  the  spectral  domain  to  solve  for 
these  unknown  coefficients,  which  in  turn  give  the  answer  to  the  unknown 
scattered  fields.  The  scattered  fields  for  all  modes  obtained  in  this 
manner  are  then  used  to  compute  the  reflection  and  transmission  coefficients 
for  each  mode,  propagating  or  attenuated.  A  tacit  assumption  made  is  that 
the  scattered  field  on  the  open-ended  side  of  the  waveguide  is  the  same  as 
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that  in  an  infinitely  long  waveguide  containing  the  beam  equalizer.  In 
other  words,  the  truncation  effects  of  the  waveguide  are  ignored  in 
this  analysis.  The  transmission  coefficients  are  used  to  weight  the 
radiation  field  due  to  each  mode  of  waveguide  and  the  superimposed 
radiation  pattern  is  computed.  The  reflection  coefficients  are  used 
to  assess  the  impedance  matching  performance.  Numerical  results  indicate 
that  the  E-  and  H-plane  principally  polarized  patterns  are  equalized 
extremely  well  over  the  entire  frequency  band  of  operation  and  that  the 
Impedance  matching  is  also  quite  satisfactory. 

II.  Analysis 

The  geometry  of  the  waveguide  with  septum  and  post  is  shown  in 
Figure  1.  Since  the  cross-section  of  the  post  is  very  small,  the  post 
is  modeled  as  a  narrow  strip  to  simplify  the  analysis.  The  incident 
field  is  propagating  in  the  z-direction  towards  the  post  as  shown 
schematically  in  Figure  1.  There  are  surface  currents  induced  on  the 
septum  and  the  post  due  to  the  incident  field.  The  scattered  fields 
radiated  by  these  induced  surface  currents  then  propagate  in  both  the 
z-direction  and  the  -z-direction,  giving  rise  to  the  transmitted  and 
the  reflected  waves,  respectively.  In  the  following  analysis,  the 
truncation  effects  of  the  waveguide  at  z»0  are  ignored,  as  though  the 
post  and  septum  were  located  in  an  infinite  guide. 

The  incident  field  in  the  waveguide  can  be  expressed  in  terms  of 
TE  and  TM  modes  in  the  usual  manner: 

TE  modes : 
nm 

K  “  Jhr>m  - r“  C0S  X^  Sin^  y)  exP  ^ 

x  nm  k  a  o  om 

c 


0  z  z  x 


LONGITUDINAL-SECTION  CROSS-SECTION 

a  =  l.3  X 
b  =  1.3  X 

c  =  0.25  X  f  =3.95  GHz 

d  IS  VARIABLE  x  =0.076  m 
t  =  0.062  X 


Figure  1.  Geometry  of  the  waveguide  with  septum  and  post. 
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H1  -  -E*  /  Z__ 
x  y  TM 


H±  -  E1  /  Z^ 
y  x  TM 


(lb) 


where  k2  -  (^L) 2  +  (S2L) 2 
c  a  o 


nm 


k2  -  k2  • 

C  » 


(a2  pe 


^TE  *  “u/Snm  ! 


S.  -  6n»/(“£> 


The  scattered  fields  can  be  expressed  in  terms  of  their  Fourier 
spectra,  which  are  in  turn  related  to  the  Fourier  spectra  of  the  induced 
surface  currents.  The  Fourier  spectra  of  the  induced  surface  currents 
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are  then  solved  for  by  the  moment  method  applied  in  the  transformed, 
i.e.,  the  spectral  domain*  Specifically,  Galerkinrs  method  is  used 
in  the  present  analysis  —  the  same  basis  functions  are  used  as  testing 
functions  in  the  moment  method.  Upon  solving  the  spectra  of  the 
induced  surface  currents,  the  scattered  fields  can  be  obtained  in 
a  straightforward  manner.  The  analytic  details  follow, 

g 

Corresponding  to  each  incident  mode,  the  scattered  fields  and 

g 

E  can  be  represented  in  the  following  form: 
z 


E®  (b/2  <  y  <  b) 
Es  (0  <  y  <  b/2) 

X 


co 

-  cos  x)  /  f  (a) 
a  '  n 


sin[y  (b-y)] 
n 


sin  Yn  y 


exp(-jaz)da. 


(2) 


"e®  (b/2  <  y  <  b) 
z 

00 

sin [y  (b-y) ] 

!e®  (0  <  y  <  b/2) 

■  sin  (~  x)  /  g  (a) 

a  4  n 

—CO 

sin  yq  y 

c  _ 

exp(-jaz)doi 


where 


1/2 


Y  -  [k2-  (ZL)2  -  a2] 

n  a 


(3) 


and  f  (a) ,  g  (a)  are  the  unknown  Fourier  spectra  to  be  determined.  The 
n  n 

s  ^ 

expression  for  E^  is  then  obtained  from  the  Maxwell's  equation  7  •  E  *  0, 

giving 


"e®  (b/2  <  y  <  bT 

00 

*  sin(^  x)  /  h  (o) 

r- 

cos(y^ (b-y)] 
n 

E®  (0  <  y  <  b/2) 

a  *  n 

mOO 

-cos  Yn  y 

exp(-jaz)da 

(4) 
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fglF"1 


with 


h  (a) 
n 


<T>  *„<«>  +  >  *„<»> 


(5) 


Now  that  we  have  E,  we  can  find  H  by  using  the  curl-of-E  Maxwell's  equation. 
However,  in  anticipation  of  relating  the  H-field  with  the  induced  surface 
currents  on  the  septum  and  the  strip,  only  the  x-  and  the  z-  components  of 
the  H-field  are  computed,  giving 


H*  (b/2  <  y  <  b) 


H*  (0  <  y  <  b/2) 


3^  sln(ir  x)  /" 


t“Y  g_(a)  +  ja  h  (a)]  cos  [y  (b-y)] 
n  n  n  n 


[y_  g„(a)  -  ja  h  (a)]  cos  y  y 
u  n  n  n 


da 
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The  induced  currents  on  the  septum  and  the  strip  can  also  be  expressed 
in  terms  of  Fourier  spectra.  For  each  incident  mode,  we  have 


J  (x,y  *  b/2,  z) 
z 


J,  (x,y  -  b/2,  z) 


sin(—  x)  /  j  (a)  exp(-jaz)  do 

3>  Z 


cos (“”  x)  /  Jx(a)  exp(-jaz)  da 


(7) 


where  j  (a)  and  j  (a)  are  Fourier  spectra  of  the  induced  surface  currents 

Z  X 

to  be  related  to  the  spectra  of  the  scattered  H-fleld.  This  relationship 
is  obtained  by  enforcing  boundary  conditions  on  the  septum  and  the  strip 
and  can  be  written  as 


~Hx  (y  '  eS  (b/2  +  e))  +  Hx  (y  "  e-*0  (b/2  •  e)) 
Hz  (y  ”  "5  (b/2  +  e))  '  S2  ^  (b/2  ’  e)) 


(8) 


where  e  is  a  positive  quantity. 

Substitution  of  (6)  and  (7)  into  (8)  leads  to  the  following  algebraic 
equations : 


Y_b 

coa  (— )  lrn 


g  (a)  -  jah  (a)] 
n  n 


j2(a) 


-2 

ja>u 

t 


Ynb 


cos  (-f-)  [y  £  (a)  +  (^)  h  (a)] 
l  an  an 


i  (o) 

X 


(9) 


I 
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Substituting  h  (a)  from  (5)  into  (9)  and  manipulating  the  resulting  equations 
n 

leads  to  the  following  matrix  equation  for  fn(a)  and  (ot>  in  terms  of  the 

transform  domain  currents  j  (a)  and  j  (o) . 

x  z 


2  2*1 

1—1 * 

(21)  k2  -  2— - 

£  (o) 

j  (o) 

a  2 

a 

n 

_ju>VYn 

“  Kb 

k2  -  a2  jo  £2-) 

8n<a) 

2cos  (  y) 

j  (a) 

X 

J 

w  -J 

(10) 


An  inversion  of  (10)  gives 


mm 

r 

.  2  2 

r  ~ 

f  (a) 
n 

-j«^) 

t2  n  ir 

j2(a) 

-j^u 

k  ~  2 
& 

gn(a) 

-  -  Y  b 

2k  t  cos  (-f-) 
n  2 

i 

CM 

1 

CM 

jx(a) 

w  - 

J 

(ID 


If  the  Fourier  transforms  of  the  currents  j  (a)  and  j  (a)  are  known,  we 

2  X 

can  obtain  fn(at)  and  gn(ct)  from  (11).  Subsequently,  (2),  (4)  and  (6)  can 

be  used  to  derive  the  scattered  fields  by  substituting  for  f  (cs)  and  g  (a) 

n  n 

in  those  expressions.  The  remaining  task  is  to  solve  for  the  transform 
currents  j  (a)  and  j  (a),  using  the  Galerkin's  method  applied  in  the 

z  X 

transform  domain. 

First  we  express  the  corresponding  space  domain  currents  in  terms  of 
a  linear  combination  of  a  set  of  suitable  basis  functions  with  unknown 


coefficients. 
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10 


v 

Qj^a)  -  ~  |  qz(z)  exp  (joz)  dz,  £  -  1,  2 . 

j  -c 


Replacing  (13)  in  (11),  one  obtains  the  expressions  for  the  transform 

domain  functions  f  (a)  and  G  (a) : 

n  n 


Fn<«> 


8n(a) 


(k2  -  i£d)  [  jQ  -  jaf  |  B^Cajj 


*  2k2Yn  cos  &£) 


jar  ^Viwj  +  <k2 


-a2)  f  i  : 


W“> 


Substituting  these  expressions  in  (2) ,  (4)  and  (6)  we  obtain  the  final 
results  for  the  scattered  fields  in  terms  of  the  unknown  coefficients 
Afs  and  Bfs.  In  order  to  find  these  unknown  coefficients,  we  enforce  the 
boundary  condition  for  the  electric  field  on  the  septum  and  the  strip, 
which  requires  the  total  tangential  component  of  the  electric  field  to  be 
zero,  i.e., 


at  y  ■  2 


Using  (1),  (2),  and  (15)  in  (16)  we  derive  the  following  equations  for  A*s 
and  Brs. 


(k2  ~  (~)  1  £  A.f*  7  tan  (Ynb  P  (at)  exp(-jaz)  da 

L  Ji-0  iJ-Yn  l2  J  1 

L  r®  y  b 

“J  ~~  £  B  I  ■£-  tan  (-?-)  Q.  (a)  exp(-jaz)  da 

Jl-1  Z  J-°>  Yn  L  % 


e  sin  (^)  exp  (+  jg  z),  for  TM  modes . 
^z  nm  z  nm 

c 


.  ^ 

2  r--  — r  h  sin(^r)  exp  (+  jg  z),  for  TE  modes. 

d  ,  z  on  z  nm 

kc 


I  «  Y  b 

j  “  E  A±  /  —  tan  <rj~)  ?±(ql )  exp  (-jaz)  da 

i»0  -«  yn 


L  -,,22  y  b 

+  £  B ^  tan  (~)Qjl(a)  exp  (-jaz)  da 

lml  —  Tn 


7  2;jk  .  ,mTi 

+  — ^ —  e  sin  (— 


T~  e  sin  (— )  exp  (+  jB  z),  for  TM  modes 
y  nm  Z  nm 


,  for  TE  modes 


Now  multiply  both  sides  of  (17)  by  the  basis  functions  p^(z)  and  q^(z)  and 
integrate  over  z  to  obtain  the  following  equations. 


2  2  2  I  *  yb  * 

(k  z  A  /  tan  (-^-)  P  (a)  P  ,  (a)  da 

a  i*0  --  Yn 


'  j  T  B*  /  T"  tan  (“T>  Va)  pi*(a)  da 
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f  -2  (~)  8 

a  rtm 

h>U 


v. 


2(f) 


k  .  mtr 

e«™  sin  T 
t  Z  nm  z 

k 

c 


k  ,  ,  mir 

-z  h  sin  -T- 

2  nm  2 

c 


pf’  «»> 

v  (6m> 


p!-  «»> 

pf  <V 


,  i’=0, for  TM  modes 


,  for  TE  modes 


i  T-  *  Ai  /  f-  tan  <Yn2)  Pi(0)  Q?’(a) 


da 


i-0  -®  Tn 


L  «  ,,2  2  Y  b 

+  ^  tan  (-y-)  Qz(a)  Q*,  (a)  da 


0 


2ik  .  inir 

—* 1 —  e  sin  — 
W4  nm  2 


QJ.  «U> 


nm 


Of  ‘V. 


,  £f«l,2,...,L  for  TM  modes 


,  i,*l>2f.*,,L  for  TE  modes 


(18) 


The  upper  quantities  within  the  square  bracket  in  (18)  are  associated 
with  positive  z  incident  modes,  and  the  lower  quantities  are  associated 
with  negative  z  incident  modes.  Both  cases  are  included,  because  we  are 
interested  in  the  transmission  coefficients  as  well  as  the  reflection 
coefficients.  The  integrals  are  evaluated  by  numerical  integration,  and 
the  resulting  system  of  linear  equations  Is  solved  as  usual  for  the 
unknowns  Afs  and  B's  by  matrix  inversion*. 
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The  choice  of  basis  functions,  i.e.,  the  p' s  and  qfs,  is  based  on 
previous  experience  and  the  consideration  of  the  behavior  of  the 
currents  at  the  edges  of  the  septum  and  the  strip.  Only  two  terms  for 
the  expansion  of  the  currents  on  the  septum  are  retained,  which  are 
believed  to  be  adequate  for  this  analysis.  The  basis  functions  and 
their  Fourier  transforms  are  shown  in  Table  1« 

Let  us  now  consider  the  case  of  interest,  namely,  one  in  which  the 
incident  fields  are  TE^  modes.  We  have  n*0,  and  the  following  relations: 


y 


2 

0 


2 

-  a 


(19) 

(20) 


When  n»0,  (18)  can  be  simplified  to  the  following 

Yob 

2  “  tan  ) 

Z  A  /  v~~K  Ma)  da 

i=»0  1  -co  Yo  1 


4  ,  * 

—  sin(-r-)  h 
mir  2  0m 


pi*’  <«W 
.Pi'  «W. 


,  i’»0,  1,  2 


2  ®  y  b 

z  B£  /  Yq  tan  Qft(a)  Q£,(a)  da  *  0,  (21) 


The  equations  for  A's  and  B's  are  uncoupled.  Therefore,  A's  and  B's  can 
be  solved  for  separately.  The  system  of  equations  for  B's  is  homogeneous. 
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Table  1.  BASIS  FUNCTIONS  IN  SPATIAL  AND  SPECTRAL  DOMAINS 


Spatial  Domain 


z+c+d+  ^  2 
pQ(z)  -  exp[-( - - - -)  ] 


Spectral  Domain 


P0(«)  "  — exp  [-ja(c+d+  j)]  exp(-(^p-)  ] 

u  4/ir  1 


px(z) 


z  +  c/2 
c/2 


) 


Pl(a)  *  t  exp  Jq(^) 


P,<*> 


z  +  c/2 
c/2 


z  +  c/2./ 

“T7r-) 


P2(a)  *  ^  exp  ^ 


^(z)  -  /l  -  0 


z  +  c/2 
c/2 


) 


Q^o) 


J  (— ) 

c  ,  ,  ca  l''  2  ; 
4  exp  (‘j  ~2)  ST 


q,(z) 


c/2 


z  +  c/2  /  ,z  +  c/2 


c/2 


-) 


J  <— ) 

n  /  \  1c  /  ,  cav  2  2  ' 
Q2(a)  -  ^  exp  (-j  — )  — - 


Jg,  J^,  and  are  Bessel  functions  of  the  first  kind  and  of  order  zero,  one,  and 
two,  correspondingly. 
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which  leads  to  the  conclusion  that  all  B's  are  zero.  Consequently,  there 
are  no  components  of  the  current  in  the  z  direction  either  in  the  septum 
or  the  strip.  In  the  following  development,  only  the  system  of  equations 
for  A's  is  investigated.  Observe  that  the  solution  of  Afs  involves  an 
inversion  of  a  3  x  3  matrix,  which  is  an  easy  task  for  the  computer.  How¬ 
ever,  the  evaluation  of  the  matrix  elements  involves  nine  complex  integrals 
to  be  numerically  integrated.  Since  there  are  singularities  in  these  integrals, 
we  must  examine  the  integrands  carefully  to  make  sure  that  the  numerical 
integration  is  applied  correctly  to  give  accurate  results  in  spite  of  the 
singularities.  Therefore,  it  is  useful  to  write  the  expressions  appearing 
in  the  integrands  of  (21)  in  an  explicit  manner  as  follows: 


P0P0  *  16ir  exP  [”2  <  4  )  ] 


P  P* 
10 


ct 


16/7” 


JQ  (y-)  exp  [-(SL)  ]  exp  [Jo(f  +  d  +  ^)  ] 


V*  J0  <f>  exp  ]  exp  H«Cf  +  *  +  !>] 


P2Po  *  J  Ji  (T>  exP  ]  exp  tja(2  +  d  + 


VS  *  *3  exp  1  exP  M«(f  +  d  +  2)] 


P1PJ  ’  16  J0  <? > 


P  p*  ■  —  j2  /£2L) 

T2  16  J1  '  2  ' 


(22) 
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Note  that  P^P*,  P2P1*  etc*  are  not  interest  here  because  the  corresponding 

integrands  are  odd  and  therefore  the  integrals  involving  them  are  identically 
zero.  Now  observe  that  the  integrand  in  (21)  contains  a  simple  pole  located 
at  a  ■  8  on  the  path  of  integration  for  0  <  a  <  «  ,  provided  that  the 
wave  number  satisfies  the  following  condition: 


(23) 


This  is  the  case  when  the  incident  field  is  the  dominant,  propagating  mode 
in  the  waveguide.  Introduction  of  some  loss  in  the  medium  clarifies  the 
position  of  the  poles  along  the  integration  path,  as  shown  in  Fig.  2.  Each 
of  the  integrals  in  (21)  can  be  written  in  the  following  form 

F(a)da  3  I  ^  F(a)da  +  [  ^  F(a)da  +  [  F(a)da  +  ir j  Res  (-g  ) 

CO  j- m  J  _G  Jo  01 


J01 


J01 


-  irj  Res  (S01) 


r  23 


01 


t 


F(a)da  +  +  F(a)da 

2801 


+  l  01  F(a)dct  + 
J0 

r-  23 


r 

hsoi 


F(o)da  +  irj  Res  (-BQ1)  -  irj  Res  (gQ1) 


01 


F(a)da  + 


B01 


[F(-26qi  'a)  +  F(a)1  da 


t 


01 


(F(ct)  +  F(23Q1  -  a)]  da 


F(a)  da  +  irj  Res  (-6^)  -  ffj  Res  (gQ1) 


01 


(24) 


J 
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In  the  above  equation,  a  bar  across  the  integral  sign  means  principal  value 
integration.  The  method  of  foldover  as  shown  makes  the  new  integrand  remain 
bounded  at  the  singularity  of  the  old  integrand;  hence,  this  new  integral 
is  easily  evaluated  by  the  numerical  method.  Also,  since  the  integrand 
goes  to  zero  rapidly  as  a  becomes  large,  the  integration  limit  <»  can  be 
replaced  by  a  large  number,  e.g.,  14  8^^.  It  should  be  noted  that  if  the 
foldover  method  was  not  used,  the  integral  could  still  be  evaluated  numer¬ 
ically,  in  some  cases,  but  the  integration  limit  00  must  be  replaced  by  a  much 
larger  number  which  requires  increased  computer  time  and  the  results  are  less 
accurate.  The  residues  given  in  (24)  are  easily  calculated  and  are  given 
by  the  general  form 

Res(±S01)  "  b^(±8m)"  Pi  (a)  V  <a)  +a  ,  i,  i'  =  0,  1,  2 

U1  01  -P01 


with  Pi  (a)  P  *  (a)  given  in  (22). 


Having  discussed  the  evaluation  of  the  matrix  elements  in  detail,  we  can 

proceed  to  solve  the  matrix  equation  (21)  to  obtain  the  A's.  Having  found 

these  Afs,  we  compute  fg(a)  from  (15)  and  then  calculate  the  scattered  fields 

s  s  s 

from  (2).  For  n-0,  both  E^  and  *  0,  and  there  is  only  E^,  which  can  be 


written  as 


K  <!  <  *  «  b> 


Es  (0  <  y  <  ^) 
x  2 


ex?  (-jaz)  da 


~j,,0llQ  I  AiPi  (a) 

2yq  cos  fy0bl  i-0 


sin  [Y()  (b  — y)  ] 


sin  yq  y 


where  yr 


-  if  k>  a 


j/a^-k^  if  k<|a| 


U 

Let  us  evaluate  (26)  at  y  =  y  , 


r  T0 

S/  b,  "jkbZ0  ?  .  ®  tan  — 


E^-f>  -  — 


l  A  /  — — r —  P.  (ot)  exp(-jctz)  da 

_ n  ^  »  nD  ^ 


1-0  -«  _0 
2 


Since  the  scattered  fields  on  both  sides  of  the  beam  equalizer  are  given  by 
(27),  we  can  compute  the  transmission  and  reflection  coefficients  by 
normalizing  these  scattered  fields  to  the  incident  field  given  in  (1). 

The  reflection  coefficient  R  can  be  expressed  as 


“*  exp(j3Q1z) 


I  A  [ 2ttJ  I  Res  (?  )  ] 
i-0  m-1,3,5, . .  « 


where  the  poles  are  given  by 


;1  *  _301 


‘2-£>  * 


r  *  i  /  C— ) 

Sn>3  2  (  b  ; 


and  the  residues  at  these  poles  are  given  by 


Resi(c1) 


4  P^-B^)  exp  (jzS01) 

»2 


4  p< (jx„)  exp  (x  z) 

Res  (c  )  =  Res  (jx)  -  - - — ^ - 5! —  ,  m  ^  3 

i  tn  l  m  2  „ 


b2(jxm> 

m 


with  x-  *  -k2 

tn  0 


PQ(jXm)  =  — ~  exp  [xm(c+d+|-)]  exp  [  (.—)  ] 
4/rr  u 


Pl(jxm)  "  t  eXp  ^  X0  (f  xm) 


P2(Jxm>  =  f  SXP  (~b  jIl  (f  V 


and  1^,  1^  are  modified  Bessel  functions  of  the  first  kind.  In  deriving  (28), 
the  integral  in  (27)  has  been  evaluated  by  the  residue  theorem  with  the 
contour  closed  in  the  upper  half  of  the  a-plane.  By  the  same  token,  the 
transmission  coefficient  T  can  also  be  obtained,  except  this  time  the  contour 
is  closed  in  the  lower  half  of  the  ot-plane. 


x  z>0 


ffexpCjS  z)  2 

-  ■  Z  A.  [-2irj  Z  Res  (c  )  ] 

4h01  i-0  1  m-1,3,5,..*  1  m 


where  the  poles  are  given  by 


r  m  8  »  j^2  ~  (JL) 

51  B01  V 


r  .  -1  /(SI)-k2 

S»3  J  r  v  b  ’  K 
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and  the  residues  at  these  poles  are  given  by 


Resi('1) 


4  P,(6q1)  exp  (-JBq1z) 


b28 


01 


.  .  ,  %  4  p^C-jXm)  exp  (-x  *) 

Res.  ('  )  *  Res,  <-Jy  )  =  - i  m  m 

i  m  i  m 


m  >  3 


with  x 


m 


b 

in 


k2 


cx_ 


P0  =  “7  exp  I_xm(c+d+  2)]  exp  l(~4S)] 

4/ir 


P1  ('V  ■  t  “f  (--r)  Xo  (f  *m> 


P2  «*P  <-x,f)  I-JIj  <fx»)l 

and  IQ,  I1  are  nodified  Bessel  functions  of  the  first  kind.  Numerical 
results  indicate  that  minimum  reflection  and  maximum  transmission  can  be 
achieved  if  the  separation  between  the  septum  and  the  strip  is  0.1 
wavelength  for  the  given  dimensions  in  Fig.  1. 

The  radiation  pattern  of  the  waveguide  feed  with  beam  equalizer 
can  now  be  computed  in  the  following  manner.  First,  the  aperture  field 
distribution  in  the  plane  containing  the  waveguide  mouth  is  estimated  by 
a  superposition  of  the  waveguide  mode  field  at  this  plane  with  each  mode 
being  weighted  by  the  corresponding  transmission  coefficients  for  the  modes. 
Since  the  transmission  coefficients  have  been  evaluated  at  a  different 
reference  plane,  it  is  necessary  to  refer  these  transmission  coefficients 
back  to  the  aperture  plane.  This  is  done  by  multiplying  the  transmission 


A 


coefficients  by  an  appropriate  correcting  factor,  which  is  a  phase  factor 
for  a  propagating  mode  and  is  an  exponential  factor  for  an  attenuated 
mode.  The  radiation  pattern  is  then  obtained  by  the  familiar  Fourier 


transform  relation  between  the  far  field  and  the  aperture  field.  For  TE 

nm 

modes,  the  far  fields  are  given  by  the  following  expressions  [23 : 

.  .  (i)1/2  il&Ltm*  (1  +  %  cos9  +  R  (1  -  %  cosQ)] 

9  e  2A'iR  lc  *  * 

nm 

*  simji)2  -  cos<ji)2]  CQ »<{>) 

a  d  nm 

1/2  2 

£  (ii)  (^rab)  sinQ  sincfr  cos 

*  G  2  A3  R 

g  g 

•  [cos©  +  +  R  (cose  --0s)]  V  (0,<J>)  (30) 

k.  k  nm 


with 

^  (9,<t>)  - 

nm 

sin(^~  sin0  cos<f>  +  — .) 

sin(“^-  sin0  sin$  +  —■) 

2  2 

(~  sin©  cos<>)  -  (~) 

(—  sin9  sin4>)  -  (^jp) 

W  W 

—  md 

.  exp  {  -j  [kR-T  sin©  (a  cos6  4*  b  sin£)  -  (n+nrhl)  ^ 

where  (R,0,<j>)  are  the  conventional  right-handed  spherical  coordinates.  The 
total  radiation  pattern  is  then  obtained  by  a  superposition  of  these 
individual  mode  patterns  with  the  appropriate  transmission  coefficients 
referenced  at  the  aperture  plane*  Computed  patterns  and  measured  results 
are  presented  in  the  next  section. 
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III.  Computed  Results 

The  performance  of  the  square  waveguide  feed  with  beam  equalizer  is 
evaluated  over  a  frequency  band  of  operation.  The  computed  radiation 
patterns  in  the  E-  and  H-planes  are  presented  in  Figure  3.  The  beam 
equalization  is  quite  satisfactory  over  the  entire  band.  The  corresponding 
aperture  field  distributions  are  shown  in  Figure  4.  Looking  at  these 
aperture  distributions,  one  can  explain  how  the  beam  equalizer  works.  It 
goes  as  follows.  The  aperture  distribution  for  the  E-plane  pattern 
is  uniform  which  gives  a  familiar  (sin  x)/x  type  pattern.  The  H-plane 
pattern  is  due  to  a  cosine-taper  type  of  aperture  distribution  if  there  is 
no  beam  equalizer  present.  Hence,  the  beamwidth  is  larger  than  that  of 
the  E-plane  pattern  for  a  square  waveguide  aperture.  However,  the  intro¬ 
duction  of  the  beam  equalizer  forces  the  H-plane  aperture  field  to  vanish 
at  the  center  of  the  aperture,  which  makes  the  field  distribution  look  more 
uniform.  Hence,  the  main  lobe  of  the  H-plane  pattern  narrows  to  achieve 
the  beam  equalization  effect  in  the  E-  and  H-planes.  The  pattern  is  mainly 
determined  by  the  septum;  the  post  is  present  for  impedance  matching.  The 
computed  reflection  coefficients  over  the  frequency  band  are  shown  in 
Table  2. 


IV.  Conclusions 

A  septum  and  an  impedance  matching  post  used  as  a  beam  equalizer  in  an 
open-ended  waveguide-feed  for  reflectors  used  in  satellite  communications 
systems  have  been  analyzed  by  using  a  spectral  domain  approach.  The 
computed  radiation  patterns  in  the  E-  and  H-planes,  as  well  as  the 
impedance  match  results,  have  been  presented  in  the  paper.  The  performance 
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Radiation  patterns  in  E-  and  H-planes  at  center  frequency  of 
a  10%  center  frequency  band. 
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Aperture  distribution  at  upper  bandedge  of  a  10%  center 
frequency  band* 
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TABLE  2 

REFLECTION  COEFFICIENTS  OVER  THE  FREQUENCY  BAND 


1 

Frequency  (GHz) 

Reflection  coefficient  r 

i 

VSWR  -  l±ill 

l-l  r  | 

3.75 

I 

0.2729/-130.250 

1.751 

3.85 

0.1879/-I34.90® 

1.463 

3.95 

0.1071/-142.02® 

1.240 

4.05 

0.0358/-167.87® 

1.074 

l  _  _ 

4.15 

0.0442/73.10° 

1.092 
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of  the  beam  equalizer  over  the  entire  band  of  operating  frequency  has 
been  evaluated.  The  results  indicate  that  the  E-  and  H-plane  principally 
polarized  patterns  are  equalized  extremely  well  over  the  entire  frequency 
band  of  operation  and  that  the  impedance  matching  is  also  satisfactory. 
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